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In any Abelian gauge theory with an action periodic in the link variables one can perform a duality transforma- 
tion not only in the partition function, but also in correlation functions including Polyakov loops. The calculation 
of expectation values in the confinement phase, like electric field strength or monopole currents in the presence 
of external charges, becomes significantly more efficient simulating the dual theory. We demonstrate this using 
the ordinary Wilson action. This approach also allows a quantitative analysis of the dual superconductor model, 
because the dual transformed U(l) theory can be regarded as limit of a dual non-compact Abelian Higgs model. 
In this way we also try to interpret the behaviour of monopole condensate and string fluctuations. Finally we 
present some applications for simulating the dual U(l) gauge theory. 



U(l) lattice gauge theory in four dimensions 
undergoes a phase transition at strong coupling 
and therefore has been widely used to investigate 
the confinement mechanism. In U(l) theory the 
condensation of magnetic monopoles is responsi- 
ble for the confinement of electric charges. The 
results of lattice simulations can be interpreted 
in terms of the dual superconductor picture of 
confinement, e.g. the validity of a dual London 
relation has been checked jl]j|]. 
On the other hand it was realized many years ago 
|U that one can perform a duality transformation 
of the path integral. In this way a new partition 
function is obtained which can be regarded as a 
certain limit of a dual non-compact Abelian Higgs 
model 0| . We will extend the duality relation 
to expectation values in the presence of external 
charges, which should help to clarify the connec- 
tion between dual transformed U(l) theory and 
the dual superconductor model. 
Our aim is to compare simulations in the dual 
theory and in ordinary U(l) gauge theory. In or- 
der to calculate expectation values in the original 
gauge theory with sufficient accuracy we use a 
very modest lattice size of 8 3 x 4. The results 
discussed below demonstrate that in the confine- 
ment phase the simulation of the dual theory is 
much more efficient, and calculations on larger 
lattices can be performed easily. 



Expectation values of physical observables O in 
the presence of a static charge pair at a distance 
d are determined by the correlation function 



{0(x)) QQ = 



(L(0)L*(d)O(x)) 
(L(0)L*(d)) 



(1) 



where L(f) is the Polyakov loop and the angle 
brackets denote the evaluation of the path in- 
tegral using the ordinary Wilson action. The 
quantities we are interested in are the electric 
field strength and the monopole currents as well 
as their squares. For the field strength we use 
the identification a 2 eF, 



sin Q^y which can be 
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shown to fulfil the Gauss law for electric charges 
Then the monopole currents can simply be 
constructed via the dual Maxwell equations. 
The results verify the dual superconductor sce- 
nario: The electric flux between a charge pair is 
squeezed into a flux tube, encircled by monopole 
currents acting like a coil. For the squared 
monopole currents, calculations on a Euclidean 
lattice show that both spatial and temporal com- 
ponents are suppressed in the flux tube. The sign 
of the spatial components, however, has to be 
changed when going back into Minkowski space 
J2J]. A profile of monopole currents and squares 
of currents at the coupling [3 — 0.96 is depicted 
in Fig. [ij to be compared to the dual theory. 
How can the obtained results be understood in 
terms of the dual superconductor picture? A very 
general starting point is the non-compact Abelian 
Higgs model. Electric charges can be introduced 
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Figure 1. Comparison of U(l) (at (3 = 0.96) and the 7 — > 00 limit of the dual Higgs model: Transverse 
profile of the monopole current (left figure) and the squared spatial (J^ 4 ) and temporal (J^ ) monopole 
currents (right figure) in the symmetry plane for charge distance d = 3a. The error bars for the results 
obtained in the dual theory have been omitted because they are much smaller than the symbols. 



by means of dual Dirac strings. The action in the 
limit of very large Higgs potential (extreme type 
II superconductor) reads 

S = 13 Gi^x^v + 2nn X!l j, v ) - 

7E^ (Ke-^^^ + h.c), (2) 

where the function G(x) determines the non- 
compact gauge field action (G(x) — > x 2 /2 in the 
continuum limit), n XlfJlV is fixed and describes 
the dual Dirac string connecting the charge pair. 
The complex Higgs field is constrained to $ T = 
e-xjp(ixx)- The ratio has the meaning of 
the squared London penetration length A in this 
model. The Higgs current J = 7/V/? (&X ~ is 
constructed in such a way that the fluxoid quan- 
tisation F + X 2 dJ = m e is valid for every single 
field configuration. The integer two-form m char- 
acterizes the physical fluxoid string which is al- 
lowed to fluctuate in simulations due to the com- 
pactness of the phase of the Higgs field. 
Let us now briefly review how to connect this 
model to the dual transformed path integral of 
U(l) gauge theory including external charges. 
Starting with the Wilson action and performing 
a Fourier expansion of the Boltzmann factor || 
one obtains the electric Gauss law as a constraint 
which is fulfilled by introducing integer link vari- 



ables *p and plaquette variables *n on the dual 
lattice. One arrives at the dual partition function 

Z= (2,) 4W I|E e-%^,103). (3) 
*p 

This corresponds to the 7 — > 00 limit of the Higgs 
model @ written in dual variables, with G(x) 
determined by the modified Boltzmann factor in 
eq. (|) and f3 dua i = 1 / '(Att 2 {3 v (1) ) . The external 
electric current becomes the boundary of a dual 
Dirac sheet *n in the dual theory. 
Finally we also perform the duality transforma- 
tion on the correlation function (|l|) to evaluate the 
expectation value of the electric field strength. It 
agrees with the electric field in the limit of the 
dual Higgs model: 

^ p _ dual \ Higgs (7^00) ,.s 

<E> QQ - <E> QQ ■ W 

Due to the validity of Maxwell equations also 
the expectation values of the monopole currents 
agree. A comparison of numerical results can be 
seen in Fig. [j]for charge distance 3a. However, an 
identity like eq. (Q) does not hold for expectation 
values of squared fields and squared monopole 
currents. The right plot in Fig. [I] shows that the 
two models behave similar, but there is no quan- 
titative agreement. 

The correspondence between U(l) gauge theory 
and dual Higgs model in the limit 7 — » 00 (this 
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Figure 2. Electric field profile in dual {7(1) at 
(3 — 0.96 between single and double charges in 
the symmetry plane for d — 3a. 



10 

8 
6 
4 
2 

0<^ 



1 1 

U -dual O 


i 


— U -dual-double + 


□ _ 


U -dual-triple □ 

□ 




□ 

+ 


+ 


□ + 




+ ° 

9 


o _ 
i 



d [a] 

Figure 3. Total energy as a function of charge 
distance for single, double and triple charges. 



means zero penetration length) allows for the in- 
terpretation that the finite extension of the flux 
tube is the result of fluctuations of Nambu's 
string. Due to the disagreement between squares, 
this picture might fail when looking at observ- 
ables like the energy density. 
Nevertheless, simulating the dual ?7(1) gauge the- 
ory gives us the possibilty to calculate any expec- 
tation value with significantly higher precision. 
This has two reasons: Firstly, the confinement 
phase is the weakly coupled phase in the dual 
theory. Secondly, it is not necessary to project 
the charge-anticharge state out of the vacuum: 
Charge pairs with arbitrary distance can be sim- 
ulated with equal accuracy, as well as multiple 
charges. As an example let us consider a pair 
of double charges at distance d — 3a. The elec- 
tric field profile in comparison to that of a singly 
charged pair is plotted in Fig. g. It can be seen 
that the flux tube widens for double charges. This 
is different from the results in deconfinement, 
where the field strength is simply doubled every- 
where. 

Finally we calculate the total energy of the elec- 
tromagnetic field as a function of charge distance 
for the Wilson action in {7(1) using the dual for- 
mulation. As mentioned above, this must not be 
confused with the energy of the dual Higgs model. 
A first calculation of this type was already per- 



formed in ref. [B using the polymer representation 
of the path integral in three-dimensional 17(1). 
Whereas in the deconfinement phase one gets the 
expected Coulomb potential proportional to the 
square of the charges, the slope in the confine- 
ment phase scales proportionally to the charge 
(see Fig. ||). These results differ from those of 
ref. H which were obtained in three-dimensional 
{7(1) gauge theory. 

We conclude that the dual transformed theory 
both helps to interpret the behaviour of expec- 
tation values and can be of great use for precise 
calculations in the confinement phase of {7(1) lat- 
tice gauge theory. Further calculations in the dual 
theory are in progress and will be published else- 
where. 
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